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Abstract 

We investigate the properties of interacting Q-balls and boson stars that sit on top of each 
other in great detail. The model that describes these solutions is essentially a (gravitating) 
two-scalar field model where both scalar fields are complex. We construct interacting Q- 
balls or boson stars with arbitrarily small charges but finite mass. We observe that in the 
interacting case - where the interaction can be either due to the potential or due to gravity - 
two types of solutions exist for equal frequencies: one for which the two scalar fields are equal, 
but also one for which the two scalar fields differ. This constitutes a symmetry breaking in 
the model. While for Q-balls asymmetric solutions have always corresponding symmetric 
solutions and are thus likely unstable to decay to symmetric solutions with lower energy, 
there exists a parameter regime for interacting boson stars, where only asymmetric solutions 
exist. We present the domain of existence for two interacting non-rotating solutions as well 
as for solutions describing the interaction between rotating and non-rotating Q-balls and 
boson stars, respectively. 



1 Introduction 

Solitons play an important role in many areas of physics. As classical solutions of non-linear field 
theories, they are localized structures with finite energy, which are globally regular. In general, 
one can distinguish between topological and non-topological solitons. While topological solitons 
PQ possess a conserved quantity, the topological charge, that stems (in most cases) from the 
spontaneous symmetry breaking of the theory, non-topological solitons [3] have a conserved 
Noether charge that results from a symmetry of the Lagrangian. The standard example of 
non-topological solitons are Q-balls [3], which are solutions of theories with self-interacting 
complex scalar fields. These objects are stationary with an explicitly time-dependent phase. 
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The conserved Noether charge Q is then related to the global phase invariance of the theory and 
is directly proportional to the frequency. Q can e.g. be interpreted as particle number [2]. 

While in standard scalar field theories, it was shown that a non-renormalisable <£ 6 -potential 
is necessary [S], super symmetric extensions of the Standard Model (SM) also possess Q-ball 
solutions [6]. In the latter case, several scalar fields interact via complicated potentials. It was 
shown that cubic interaction terms that result from Yukawa couplings in the superpotential 
and supersymmetry breaking terms lead to the existence of Q-balls with non- vanishing baryon 
or lepton number or electric charge. These supersymmetric Q-balls have been considered as 
possible candidates for baryonic dark matter [7] and their astrophysical implications have been 
discussed j8j. In [21], these objects have been constructed numerically using the exact form of 
the supersymmetric potential. 

Q-ball solutions in 3 + 1 dimensions have been studied in detail in [5j [TU1 [TT] . It was realized 
that next to non-spinning Q-balls, which are spherically symmetric, spinning solutions exist. 
These are axially symmetric with energy density of toroidal shape and angular momentum 
J = kQ, where Q is the Noether charge of the solution and k G Z corresponds to the winding 
around the z-axis. Approximated solutions of the non-linear partial differential equations were 
constructed in [5] by means of a truncated series in the spherical harmonics to describe the 
angular part of the solutions. The full partial differential equation was solved numerically in 
\10\ [TTJ [T2]. It was also realized in [S] that in each fc-sector, parity-even (P = +1) and parity- 
odd (P = —1) solutions exist. Parity-even and parity-odd refers to the fact that the solution is 
symmetric and anti-symmetric, respectively with respect to a reflection through the x-y-plane, 
i.e. under 6 — > ir — 8. 

These two types of solutions are closely related to the fact that the angular part of the 
solutions constructed in [5j \W\ [TT] is connected to the spherical harmonic Y®(9, (p) for the 
spherically symmetric Q-ball, to the spherical harmonic Y^(Q, (p) for the spinning parity even 
(P = +1) solution and to the spherical harmonic Y^O, (p) for the parity odd (P = —1) solution, 
respectively. Radially excited solutions of the spherically symmetric, non-spinning solution were 
also obtained. These solutions are still spherically symmetric but the scalar field develops one 
or several nodes for r €]0, oof. In relation to the apparent connection of the angular part of the 
known solutions to the spherical harmonics, "^-angular excitations" of the Q-balls corresponding 
to the spherical harmonics Y^{0, ip), —l<k<l have been constructed explicitely for some values 
of k and I in [12]. These excited solutions could play a role in the formation of Q-balls in the 
early universe since it is believed that Q-balls forming due to condensate fragmentation at the 
end of inflation first appear in an excited state and only then settle down to the ground state [13] . 
The fact that these newly formed Q-balls are excited, i.e. in general not spherically symmetric 
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could, on the other hand, be a source of gravitational waves |14j . 

Interactions of well-separated Q-balls in (1 + l)-dimensions have been studied in [15] and it 
was shown that Q-balls with equal frequencies can attract when being in-phase or repell when 
being exactly out-of-phase. 

The interaction of two Q-balls in (3 + l)-dimensions that sit on top of each other has been 
studied in |12| . It was found that the lower bound on the frequencies cjj, i = 1,2 is increasing 
for increasing interaction coupling. Explicit examples of a rotating Q-ball interacting with a 
non-rotating Q-ball have been presented. 

Complex scalar field models coupled to gravity possess so-called "boson star" solutions [TBI 
[TTl HEl UHl [20l EJJ . In [TUl[TTl[22] boson stars have been considered that have flat space-time limits 
in the form of Q-balls. These boson stars are hence self-gravitating Q-balls. The interaction of 
boson stars has been studied in [22] and it was found that ergoregions can appear when a non- 
rotating boson star interacts with a rotating and parity even boson star signaling an instability 
of the solution. Recently, charged Q-balls and boson stars in scalar electrodynamics have also 
been considered [23]. 

In this paper, we study interacting boson stars and Q-balls that sit on top of each other. 
While in [121 E2] we were mainly interested in the different types of solutions existing in the 
model, we present here an analysis of the dependence of charges and masses of the solutions on 
the parameters of the model. During this analysis, we have found that the pattern of solutions 
is much richer than in the non-interacting case, especially we observe new (and unexpected) 
branches specific to the system of two scalar fields. Fixing the different coupling constants ap- 
pearing in the Lagrangian in such a way that it is symmetric under the exchange of the two scalar 
fields, we impose an extra Z2 symmetry. Then, it turn out that the new branches correspond to 
solutions which break this symmetry. In other words, the Z2 symmetry is spontaneously broken 
by the new solutions. 

The paper is organized as follows: in Section 2, we give the model, Ansatz and boundary 
conditions. In Section 3 and 4, we discuss our results for non-rotating and rotating solutions, 
respectively. Section 5 contains our conclusions. 



In the following, we study a scalar field model coupled minimally to gravity in 3 + 1 dimensions 
describing two interacting boson stars. The action S reads: 



2 The model 
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where R is the Ricci scalar, G denotes Newton's constant and C m denotes the matter Lagrangian: 



~d^ x d^\ - \d^ 2 d^* 2 - y($i,$ 2 ) (2.) 



where both $1 and $2 are complex scalar fields and we choose as signature of the metric 
( — h ++). The potential reads: 

2 

$2) = N*i| 6 - ^l*i| 4 + MH 2 ) + 7l^i| 2 |^2| 2 (3) 
i=l 

where k%, Pi, Aj, i = 1,2 are the standard potential parameters for each boson star, while 7 
denotes the interaction parameter. The masses of the two bosonic scalar fields are then given 
by (m%) 2 = X h z = l,2. 

Along with |10|. [TT1 [TZ] , we choose in the following 

Ki = l , pi = 2 , A; = 1.1 , i = l,2 . (4) 

This particular choice of parameters leads to an extra Z 2 symmetry of the Lagrangian: $1 <-> <3>2- 
In [5] it was argued that a $ 6 -potential is necessary in order to have classical Q-ball solutions. 
This is still necessary for the model we have defined here, since we want $1 = and $2 = to 
be a local minimum of the potential. A pure $ 4 -potential which is bounded from below wouldn't 
fulfill these criteria. 

The matter Lagrangian C m ([2]) is invariant under the two independent global U(l) transfor- 
mations 

$1 -► $ie ixi , $ 2 -» <&2e iX2 ■ (5) 

As such the total conserved Noether current j? tat \, fJ> = 0, 1,2,3, associated to these symmetries 
is just the sum of the two individually conserved currents and j% with 

tftot) =j?+j$ = -i(nd^i-$id»$t)-i(^d^ 2 -$ 2 d^* 2 ) . (6) 

with ;/1 = 0, # ;M = and ;/i = 0. 

The total Noether charge Q(tot) °f the system is then the sum of the two individual Noether 
charges Q\ and Q 2 : 

Qitot) =Ql + Q2 = ~J f^X - J j$d 3 X (7) 

Finally, the energy-momentum tensor reads: 

2 

T^ u = {d^id^* + d u ^d^l) - g^C (8) 
i=i 
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2.1 Ansatz and Equations 

For the metric the Ansatz in Lewis-Papapetrou form reads |10| : 

ds 2 = -fdt 2 + L [g{dr 2 + r 2 d0 2 ) + r 2 sin 2 6 {dip + ™^) 2 ) (9) 

where the metric functions /, I, g and m are functions of r and 9 only. For the scalar fields, the 
Ansatz reads: 

$i(t, r, 9, tp) = e^ t+ik ^Ur, 0) , i = 1, 2 (10) 

where the and the ki are constants. Since we require $i(</?) = &i{<p> + 2-7r), i = 1,2, we have 
ki £ Z. The mass M and total angular momentum J of the solution can be read off from the 
asymptotic behaviour of the metric functions [10]: 

M = — — lim r 2 d r f , J = —— lim r 2 m . (11) 

2G r^oo 2G r^oo 

The total angular momentum J = J\ + J2 and the Noether charges Qi and Q2 of the two boson 
stars are related by J = k\Q\ + A/2<32- Boson stars with fcj = have thus vanishing angular 
momentum. Equally, interacting boson stars with fci = — k2 and Qi = Q2 have vanishing angular 
momentum. 

The coupled system of partial differential equations is then given by the Einstein equations 

= SirGT^ (12) 

with T^y given by ([8]) and the Klein-Gordon equations 

dV \ 

D+ d\^)^ = ° ' i = 1 ' 2 ' (13) 
The explicit expressions for the equations can be found in the Appendix. 

So far, the scale of the scalar fields is not yet fixed and the functions (fri, 4>2 are dimensionful. 
In order to study the equations, it is convenient to rescale these function according to 

fa^&r) , » = 1,2. (14) 

where rj has the dimension of an energy and its scale depends on the chosen phenomenological 
model [SI [TJ El The potential parameters are rescaled as 7 — ► 7/f? 2 , /3% — > fti/rj 2 and 

/tj > K{jT\ . 

We can then introduce the dimensionless quantity 

a = 8ttG V 2 = ^ (15) 



2 _ 87Tff 
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which measures the ratio between the energy scale of the scalar field and the Planck mass. 

As stated in [24J, in supersymmetric extensions of the standard model, the value of r/ is 
between a few TeV and a few hundred TeV. The value of a is hence very small. However, in 
order to be able to compare our results with the case of a single Q-ball or boson star, we follow 
the discussion in |10t [TT] and also study larger values of a throughout this paper. Furthermore, 
larger values of a also accommodate possible other energy scales. 

2.2 Boundary conditions 

We require the solutions to be regular at the origin. The appropriate boundary conditions read: 
d r f\ r=0 = , <9 r Z| r=0 = , g\ r=0 = 1 , m\ r=0 = , <^| r=0 = , £ = 1,2. (16) 

for solutions with k{ ^ 0, while for ki = solutions, we have d r (fii\ r= o = 0, £ = 1,2. The 
boundary conditions at infinity result from the requirement of asymptotic flatness and finite 
energy solutions: 

/Ir— >oo — 1 j l\r^oo — 1 > >oo — 1 > T7l| r _ >DO — , </>i|r— »oo — , £ — 1,2. (17) 

For 9 = the regularity of the solutions on the z-axis requires: 

def\e=o = , d e l\e=o = , g\ 6=0 = 1 , d e m\ e=0 = , c/>i\e=o = , £ = 1,2, (18) 

for ki ^ solutions, while for ki = solutions, we have de<f)i\ r= Q = 0, i = 1, 2. 
The conditions at 9 = ir/2 are either given by 

def\e=7r/2 = , d e l\ e=n / 2 = , d e g\ g=n / 2 = , d e uj\ 0=n/2 = , d e 4>i\ g=7T/2 = , £ = 1,2 

(19) 

for even parity solutions, while for odd parity solutions the conditions for the scalar field functions 
read: (j)i\ e=n / 2 = 0, £ = 1,2. 

3 Non-rotating solutions 

The solutions have vanishing angular momentum for k% = k 2 = 0. In this case, the system of 
differential equations reduces to a system of coupled ordinary differential equations and g = 1, 
m = and the remaining functions are functions of the radial variable r only. 

We have solved the corresponding ordinary differential equations (ODEs) using the ODE 
solver COLSYS [25]. 
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3.1 a = : Interacting Q-balls 

In the flat space-time limit, i.e. for a = the metric functions / = I = 1 and the system 
describes two interacting, non-rotating Q-balls. These Q-balls are interacting only if 7 7^ via 
the potential interaction. For single Q-balls, it has been observed [5J fTUl [TT] that the solutions 
exist only on a finite interval of the frequency uj. It was realized in [12] that this is also true for 
interacting Q-balls. In the limit of two non-rotating Q-balls the upper bounds on u\ and LV2 are 

u\ < u\ max = Ai , lo\ < u\ max = A 2 (20) 
To find the lower bound, we introduce a polar decomposition of <p\ and 4>2- 

(j) 1 = pcosx , <A2 = psinx (21) 
such that the lower bound reads |12|: 



tof cos 2 x + sin 2 x > (w 2 cos 2 x + u\ sin2 X) 



nun 



. > 1 (/3icos 4 x + /32sin 4 x-7cos 2 xsin 2 x) 

Aicos x + A 2 sm x~ ~, » ; — (22) 

4 m cos D x + k-2 sin x 
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We have studied the dependence of the charges and the mass of the solution in dependence 
on the potential parameter 7. We have also investigated how the mass of the solution evolves in 
comparison to the mass of Q\ + Q2 individual bosons with masses tub 1 and mg^, respectively. 

The pattern of solutions is involved when studying the dependence of the solutions on u±, 
ix>2 and 7. Thus, we first discuss the case u\ = 0J2, where we expect in analogy to the non- 
interacting case with 7 = that we should always find = <p2, i-e. Qi = Q2- However, already 
in this limit, the pattern of solutions turns out to be richer than expected. We come back to 
this observation at the end of this subsection. First, let us discuss the effect of 7 on solutions 
with Qi = Q 2 . 



3.1.1 uii = 0J2 

We first discuss the case (pi = <f> 2 , i.e. x = 7r/4. Then the bounds on the frequencies are given 
by U2DD and (see (l22l): 

K + ^)min = \+l-\l 2 <"l+"l (23) 

For lo\ = UJ2 tending to the upper bound (|20|) the charges Q\ = Q2 and the mass M diverge. 
This is clearly seen in Figs Q] and Gl where we plot the dependence of the charges Qi = Q2 and 
the mass M of the solution on u\ = 0J2 for different values of 7. Apparently the charges and the 
mass diverge for u>i — > \/Ai ~ 1.0488. 
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Figure 1: The dependence of the charges Q\ = Q2 of two interacting Q-balls on u\ = L02 is 
shown for different values of 7. The small numbers indicate the value of 7. 

For ui = LJ2 tending to the lower bound (|23p . we observe that the limiting behaviour depends 
crucially on the value of 7. To understand this, we observe that the lower bound (|23[) becomes 
negative for 7 < 4 - t/88/5 « -0.195 and 7 > 4 + ^/88/5 « 8.195. This means that for the 
corresponding values of 7 we can lower u\ = 0J2 down to zero. For lo\ = L02 = 0, the effective 
potential 

V eff = \ (uUi+u^l) - Vifhifa) (24) 

is nowhere positive and the solutions are unphysical. Nevertheless, let us stress that we can 
have solutions with arbitrary small charges (and frequencies). Note that this feature is not a 
consequence of the interaction but is also present in the single Q-ball case for an appropriate 
choice of (3 different from that used in [5j flO|. [TT] . 

It is apparent in FigsH]and[2l where we show the charges Q\ = Q2 and the mass M of the 
Q-ball solutions for negative values of 7, that solutions exist down to u\ = u>2 = 0. The limiting 
solutions have vanishing charges, but non- vanishing mass. 

We also show the binding energy M — (Qim^i + Q2fUB,2) = M — (Q\ + Q2)m,B of the 
solutions in Fig|3l This quantity compares the mass of the interacting Q-balls M with the mass 
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Figure 2: The dependence of the mass M of two interacting Q-balls on uj\ = loi is shown for 
different values of 7. The small numbers indicate the value of 7. 

of Qi + Q2 bosons with respective masses tub,! = tub and itib,2 = m B and indicates whether 
the Q-balls are stable. For negative and positive binding energy, we expect the solutions to be 
stable and unstable, respectively. We observe that for 7 > 0, the solutions are stable for nearly 
all values of the frequencies uji = L02, apart from values close to the maximal frequency. 

This changes for negative values of 7. For 7 < —2, the binding energy is positive for all 
values of lo\ = u>2 indicating an instability of the solution. For 7 = — 1, the solutions are stable 
in the interval lo\ = U02 € [0.15 : 0.7], but unstable for all other values of the frequency. 

For sufficiently strong interaction between the Q-balls and uj\ = L02, we observe a new 
phenomenon. For fixed values of the interaction parameter 7 and the frequency uj\ = UJ2, we 
find that two types of solution exist. One solution has = 4>2 ("the symmetric solution" in 
the following, see discussion above), while the second solution has 4>i 7^ <p2 ("the asymmetric 
solution" in the following). This is demonstrated in FigU]for 7 = —0.5 and lo\ = U2 = 0.6 where 
we give the profiles of <j)\ and 4>2- For the symmetric case, (pi = 4>2, while for the asymmetric 
case (pi 7^ (p2- Note that in this case, we present a solution for which (/>i(0) < </>2(0), but that 
due to the symmetry of the equations a second solution with (pi and 4>2 interchanged exists. In 
this sense, there is - in fact - not only one asymmetric solution for a fixed oji = 0J2, but two. 
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Figure 3: The dependence of the binding energy M — {Q\ + Q2) m B of two interacting Q-balls 
on u)\ = L02 is shown for different values of 7. The small numbers indicate the value of 7. 



Moreover, it is obvious that 4>i/<f>2 is not simply a constant. We observe that the asymmetric 
solutions have much higher energy than the symmetric ones, e.g. for the case plotted in FigUl 
we have E syrnrne t r i c ~ 146 for the energy of the symmetric solution, while the energy of the 
asymmetric solution is Easymmetric ra 288. For fixed 7, we observe that the asymmetric solutions 
exist only up to a maximal value of uj\ = u>2 = uj, which we call uj cr in the following. At uj cr , the 
asymmetric solutions join the branch of symmetric solutions. We find that uj cr depends crucially 
on the choice of 7. For 7 > —0.195, asymmetric solutions do not exist at all, while the critical 
value of uj increases with decreasing (and negative) 7, e.g. we find uj cr (^ = —0.5) » 0.74 while 
w cr(7 — —1-0) ~ 0.82. This means that the stronger the repulsion between the Q-balls, the 
higher the frequency of the asymmetric solution can be. 

The domain of existence of Q-ball solutions with lo\ = U2 = uj is shown in Figj5] (blue 
curves). The curves for uj m i n and uj max are given by (ffi)]) and ([^]) . while the curve ui cr has 
been determined numerically. As is obvious from this plot, asymmetric solutions exist only for 
7 < -0.195. 

Let us emphasize that whenever asymmetric solutions exist, a corresponding symmetric 
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Figure 4: The profiles of the functions fa an d fa are shown for symmetric Q-ball solutions and 
asymmetric Q-ball solutions, respectively. Here u)\ = UJ2 = 0.6 and 7 = —0.5. 

solution with lower energy is also present. In a concrete physical setting, we would thus expect 
the asymmetric solutions to be unstable with respect to a decay into the symmetric solutions. 

3.1.2 uj\ 7^ U2 

In this subsection, we would like to discuss the pattern of solutions describing the repulsive 
interaction of two Q-balls with non-equal frequencies. In the case 7 = 0, the solutions exist in a 
rectangle of the Wi-u^-plane which consists of the direct product [(jJ\ im im <^i,max] x [ w 2,min! w 2,m(a] • 
This is represented by the green rectangle in Fig. [6j For 7 < (we set 7 = —1) the pattern of 
solutions changes non-trivially. Since the parameters of the potential are chosen in such a way 
that the equations are symmetric under fa <-> fa, the domain should be symmetric under the 
reflection uj\ <-> 002 ■ Accordingly the discussion of the domain is easier by using the parameters 
S = (u\ + W2)/2 and A = {uj2 — 0J\)/2. We observe that up to three branches of solutions exist 
depending on the choice of S and A. It turns out that these various branches interconnect the 
symmetric and asymmetric solutions present in the case oj\ =102- This is demonstrated in FigEJ 
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Figure 5: The domain of existence for Q-balls with u\ = 002 = w is shown in the u-7-plane for 
a = (blue) and a = 0.1 (red), respectively. 

where we give the domain of existence of two Q-balls for 7 = — 1 in the wi-u^-plane. The labels 
1, 2, 3 in the figure refer to the number of solutions, while the labels a, b, c, d are discussed in 
the following. It turns out that the domain of existence can be separated into four subdomains. 
To illustrate this, we plot 0i(O) and 02 (0) ( see FigJT]) as well as the mass M of the solutions (see 
FigJS]) for four different fixed values of E in dependence on A. Note that the domain of existence 
is given in the ui-u^-plane, but that we discuss the solutions using A and E. Constant values 
of E correspond to diagonals that fulfill uji = 2E — lo\ in the cJi-tJ2-plane. These diagonals can 
be easily followed in FigEJ 

1. For 0.96 < E < 

^max the symmetric solution <fii — (f>2 gets deformed and exists for 
-E/2 + uj max > A > E/2 — oj max . On the boundaries of this interval, the charges and 
mass of the solutions diverge (like in the non- interacting case). 

2. For 0.82 < E < 0.96 (see E = 0.9 in Fig.sElED interacting solutions exist for A b < A < A a 
(with A b = -A a ). At A = A a f> one of the two scalar fields vanishes identically and the 
limiting solution does no longer represent interacting Q-balls, but only single Q-balls. The 
values of are represented by the red lines labelled a and b respectively. 

3. For 0.53 < E < 0.82 (see E = 0.8 and E = 0.6 in Fig.sEl E]), a new phenomenon occurs. 
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When deforming the solution = <f>2 by means of A > 0, we produce a branch of solutions 
which terminates at some critical value of A, say A = A c > 0. This new critical value 
is represented by the black line and by the label "c" . A second branch of solutions then 
exist on A E [A c , A&]. On this second branch, the mass of the solutions is higher than 
on the main branch. The second branch extends in particular into the region ui = L02 
but the corresponding solution has §2 7^ <p\. This branch of solutions clearly connects the 
symmetric with the asymmetric solution for lo\ = u>2- Deformation of the solution 0i = (p2 
by means of A < is also possible and leads of course to the occurrence of a third set of 
solutions for A G [A a ,A^]. 

4. Finally, for £ < 0.53 (see £ = 0.4 in Fig.s [5]) the symmetric solutions at A = get 
deformed for A > up to a maximal value A c (limiting curve labelled "c" in Figj6|). 
Starting from this value of A a second branch of solutions exists for A 6 [A c , A a ], where 
A a denotes the value of A on the curves labelled "a". For A < 0, the pattern is similar, 
but the deformed symmetric solutions exist down to A^ and from there a second branch 
of solutions exist up to A&. 

3.2 a/0 : Boson stars 

For a 7^ the solutions are self-gravitating Q-balls, so-called boson stars. In this case the 
solutions interact both via the potential interaction as well as through gravity. 

3.2.1 uj\ = u>2 

First, we have fixed a = 0.1 and have studied the dependence of the charges and the mass on 
the interaction parameter 7 for the symmetric case (pi = <p2 (for the case 7^ <p2 see the end 
of this section). Our results are shown in Figs EE Hniand[TTJ Like in the flat space-time limit, 
the solutions exist on a finite domain of the frequency. As for single boson stars [101 lllj . the 
charges Q\ = Q2 and the mass M tend to zero for u\ — > u>i :Tnax (u>2 — ► ^2, max)- The maximal 
value of loi = u>2 seems to be practically independent of a and is equal to the flat space-time 
value. This is shown in Figs£J] and [TUl where we give the dependence of the charges Q\ = Q2 
and the mass M, respectively, on the frequency oj\ = 002- In this limit the value of the metric 
function /(r) at the origin, /(0), tends to one and since f(oo) = 1, we find f(r) = 1. This is 
plausible since in this limit, the mass tends to zero and there is thus no energy- momentum to 
curve the space-time. 

For u)\ = u>2 tending to the minimal value we observe an inspiralling of the charges as well 
as of the mass. This is very similar to the case of single boson stars |10[ 111], At the same time, 
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the value of /(0) tends to zero as can be seen in Figllll However, we find that this pattern of 
inspiralling exists only as long as the minimal value of the frequency is non- vanishing. When 
the solutions exist down to u\ = lo<i = 0, the pattern changes. For 7 = —0.5 and 7 = —1 we 
observe that the charges Q\ = Q2 tend to zero, but that the mass is non-vanishing. The limiting 
solution is thus a non-trivial, static solution. Moreover, this solution possesses non-trivial metric 
functions as is obvious from the fact that /(0) 7^ 1 (see Fig Jlip . This limiting solution is thus a 
self-gravitating static scalar field solution. 

The value of 7 for which solutions exist down to oj\ = 002 = 0, 70, depends on a. In the flat 
space-time limit, this value of 7 is given analytically by (I23p and we have 70 (a = 0) ~ —0.195. We 
find that this value of 70 decreases for increasing a. For a = 0.1, we have 7o(a = 0.1) ~ —0.49. 
This is shown in Figl5] (red curves). Interacting boson stars exist for all values of lo\ = u>2 above 
the curves. Note that the a = curve is the analytic curve given in (j23l) . 
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Figure 7: The value of </>i(0), respectively 02 (0) is shown in dependence on A for different values 
of X. The labels "1" and "2" refer to the components and 02- 
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Figure 8: The value of the mass M is shown in dependence on A for different values of X. 

Again, we observe that next to symmetric solutions, asymmetric solutions exist. A typical 
asymmetric solution is shown in FigfT2l for a = 0.1, 7 = —0.5 and u\ = UJ2 = 0.6. For 
comparison, we also plot the corresponding symmetric solution. The profiles of the scalar field 
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Figure 9: The dependence of the charges Q\ = Q2 of two interacting boson stars on u\ = lo<i is 
shown for different values of 7 and a = 0.1. The small numbers indicate the value of 7. Note 
that we plot Qi/3 = Q2/3 for 7 = —0.5. 

functions <p\ and <p2 look very similar to what has been observed in the flat space-time limit, 
however, the boson stars have smaller radii as the corresponding Q-balls for the same values of 
the parameters. The metric functions of the symmetric solution differ only little from those of the 
asymmetric solution. We observe that both /(0) and 1(0) are slightly higher for the asymmetric 
solution. That the flat space-time phenomenon persists when gravity is added is not that 
surprising. However, we have also studied this phenomenon for vanishing potential interaction 
7 = and have chosen a = 0.1 letting the boson stars with uj\ = L02 = oj interact only via 
gravity. We find that symmetric solutions exist on the interval [w^ in : w„ a3 .] » [0.456 : 1.041], 
while asymmetric solutions exist on the interval [u^f n : oj^f ax ] [0.355 : 0.462]. Thus, there is 
a small interval in which both symmetric as well as asymmetric solutions exist. In this interval 
the mass of the symmetric solution is much smaller than that of the asymmetric one. We would 
thus expect the asymmetric solutions to be unstable to decay to the symmetric ones. However, 
for 0.355 < uj < 0.456 only asymmetric solutions exist. This is a completely new phenomenon 
as compared to the flat space-time limit. 
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Figure 10: The dependence of the mass M of two interacting boson stars on lo\ = UJ2 is shown 
for different values of 7 for a = 0.1. The small numbers indicate the value of 7. Note that we 
plot M/2 for 7 = -0.5. 

When studying the binding energy of the solutions, we observe that the curves look quali- 
tatively similar to the ones shown in FigEl this is why we don't present an extra figure here. 
Quantitatively, boson stars are stronger bound than their corresponding flat space-time coun- 
terparts, which is - of course - due to the attractive nature of the gravitational interaction. 

3.2.2 oji / ll> 2 

Here, we concentrate on the influence of the gravitational interaction on the domain of existence 
of non-rotating boson stars. Let us recall that the solutions exist in a rectangle of the wi-u^- 
plane which is described by the interval u^i t max] x [^2,min, ^>2,max]- The domain of Q-balls 

in the a = limit is shown in Fig. [T3] (green lines). 

For 7 = 0, q > (we have chosen a = 0.1) the pattern changes considerably. Again, it is 
useful to discuss the solutions by using £ = (loi + u>2)/2 and A = (u>2 — U\)/2. Constant Ss 
correspond again to diagonals in the a;i-a;2-plane. In the A = limit, the symmetric solutions 
4>i = <p2 exist for S G [0.456, 1.041] and can de deformed into solutions with A ^ 0. For A — > A a 
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Figure 11: The dependence of the value of the metric function / at the origin, /(0) is shown for 
interacting boson stars for different values of 7 and a = 0.1. 

(resp. A — ► — A a ) the field <$>\ (resp. fa) vanishes identically and the limiting solution represents 
a single boson star. This is demonstrated in Fig.s I14|15l for £ = 0.6. In Figll4|. it is apparent 
that for A — > A a ~ 0.55, the charge Q\ tends to zero (see also Fig ll51 where 0i(O) — > in this 
limit). For A — > —A a ~ —0.55 we find Q2 — > 0. For Qj — ► 0, i = 1, 2, a single boson star 
remains that has finite mass and charge. 

The corresponding domain in the wi-u^-plane and the critical values ±A a are labelled by 
the symbol "S" and the lines "a" in Fig. \13\ respectively. 

For £ 6 [0.355, 0.465] and A = 0, we observe that a pair of new solutions exists bifurcating 
from the branch of symmetric solutions at £ ~ 0.465. These solutions have <pi 7^ 4>2 an d 
correspond to asymmetric interacting boson stars. One of the solutions has </>i(0) > ^2(0), while 
the other has 0i(O) < </>2(0). Note that the two solutions are, of course, related to each other 
by the exchange <p\ <-> <p2- 

Deforming the asymmetric solutions for A > 0, it turns out that the solutions which have 
(/>i(0) < 02(0) at A = will have a vanishing scalar field (j)i(r) = at Af, (denoted by the red line 
"b" in FiglED. This is demonstrated in Fig.s H US] for S = 0.4, where Qi = and ^i(0) = 

18 




Figure 12: The profiles of the scalar field functions (f>±, 4>2 and of the metric functions / and I are 
shown for symmetric boson stars and asymmetric boson stars, respectively. Here oj\ = UJ2 = 0.6, 
7 = —0.5 and a = 0.1. 



(and hence also </>i(r) = 0) at A& ~ 0.06, while the mass M, the charge Q2, </>2(0) and /(0) stay 
finite in this limit. On the other hand, if we deform the solution which has 4>i(0) < 4>j(0) at 
A = for A < 0, we find that /(0) tends to zero (after a possible inspiralling) at A c (denoted 
by the blue line "c" in Fig |13p . This phenomenon is demonstrated in Fig|15l for £ = 0.4 and 
4>i (0) < 2 (O) at A = 0. Clearly for A « -0.03, we have /(0) « 0. 

Note that if we would choose the asymmetric solution with 0i(O) > 02(0) at A = the 
pattern would be the same, but with A^ — ► — and A c — ► — A c . 

In the small domain between the curves "b" and "c" only solutions which are deformations 
of the asymmetric solution at A = (denoted by "A") exist. This is completely different from 
the flat space-time case, where asymmetric solutions exist only when also symmetric solutions 
are present for u% = u%. This means that for appropriate choices of the frequencies, boson stars 
could exist that have the same frequencies but different charges. These asymmetric solutions 
would then be the generic solutions in this parameter range. 
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Figure 13: The domain of existence of interacting boson stars in shown in the wi-u;2-plane for 
a = 0.1. 

4 Rotating solutions 

As pointed out already in previous sections, the equations for the fields 4>i and <j)% totally decouple 
in the limit 7 = 0, a = and solutions of the different types available in the case of a single Q- 
ball can arbitrarily be superposed. Solutions of the decoupled system can naturally be labelled 
according to the "quantum numbers" characterizing each of the two scalar fields, say /q, l{ for 
i = 1,2. Any such configuration then gets deformed by the direct interaction if 7 7^ and/or 
by gravity if a 7^ 0. The domain of existence in the wi-u;2-plane is rectangular when a = 7 = 
and we observe this domain to be deformed when 7 and/or a are non- vanishing. The study of 
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Figure 14: The value of the mass M and the charges Q\ and Q2 are shown in dependence on 
A = (oji — u)2j/2 for two different choices of £ = (u\ + u)2)/2. Here a = 0.1 and 7 = 0. 

the evolution of this domain turns out to be very involved. In this section, we present some of 
its qualitative features for particular cases. More precisely, we study rotating Q-balls (or boson 
stars) interacting with non-rotating Q-balls (or boson stars). We have chosen the fundamental 
non-rotating solution characterized by k\ = 0, l\ = to interact with the rotating solution with 
k 2 = l,h = I- 

We have solved the system of partial differential equations (fT2j) and (fT3|) subject to the 
appropriate boundary. This has been done using the Partial differential equation solver FIDISOL 
|26j . We have mapped the infinite interval of the r coordinate [0 : 00] to the finite compact 
interval [0 : 1] using the new coordinate z := r/(r + 1). We have typically used grid sizes of 150 
points in r-direction and 70 points in 6 direction. The solutions presented here have relative 
errors of 10 -3 or smaller. 
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Figure 15: The value of the scalar field functions and of the metric function / at the origin 
4>i (0), (^2(0), f(0) are shown in dependence on A = (lo\ — u>2)/2 for two different choices of 
£ = (cji +CJ2V2. Here a = 0.1 and 7 = 0. 

4.1 a = 0: Interacting Q-balls 

For 7 = 0, the domain of existence of the fci = 0, l\ = and k% = 1> h = 1 solutions consists of 
the rectangle 

£ [^l.min.) Wi imax ] , W2 £ [<JJ2,min, ^2,max] ■ (25) 

This is indicated by the black line in Fig. [T6j The boundary of the rectangle is given by [10] : 



'1.1 



^>l,min — yl/10 ~ W2,r, 



(26) 



In fact, on basis of numerical results, it has been conjectured in [10J that u>2 
but this equality could not be proved analytically. For both limits uji — ► ca 



i,max j 



cjj — > uJi tm i n , i = 1,2 the corresponding charges diverge. When cjj 
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Figure 16: The domains of existence of non-rotating Q-balls (boson stars) interacting with 
rotating Q-balh (boson stars) are shown in the wi-u;2-plane for different values of 7 and a. 



fields <pi spread completely over the full interval of r, while for u>i — > uii >m i n , i = 1,2, they become 
step- like. 

We now discuss qualitatively the deformation of the domain in the case of interacting Q- 
balls, i.e. for 7 7^ 0. We have analyzed in detail the cases 7 = —0.5 and 7 = — 1 which capture 
the main features. The form of the domain is given in Fig. [16] for 7 = — 1 by the line with green 
triangles. 
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Figure 17: The scalar field functions fa and fa for a system of a non-rotating Q-ball with 
k\ = l\ = interacting with a rotating Q-ball with k^ = I2 = 1 and 7 = —0.5. The plot on the 
left and right corresponds to u\ = 0.5, L02 = 0.85 and ui\ = 0.5, 002 = 0.525, respectively. 

4.1.1 Fixed oj\ 

For fixed u\ the solutions exist on a finite interval of o^- The limits of this interval depend only 
slightly on ui, but strongly on 7. E.g. for uj\ = 0.5 we find 

uj 2 G [0.52,0.85] for 7 = -0.5 , u; 2 G [0.75,0.95] for 7 = -1.0 

i.e. decreasing 7, the interval in for which solutions exist decreases and gets shifted to higher 
oj2- For UJ2 close to 

^2 mm 1 the profile of <pi hardly deviates from the spherically symmetric 
solution (i.e. fa has its maximum at the origin). The profiles of the two scalar functions are 
shown in Fig[T7] (left) for u)\ = 0.5, L02 = 0.85 for 7 = —0.5. The solution stops for a reason which 
remains unclear, it is very likely that a second branch of solution backbending from the first 
branch exists. Should it exist, it would constitute the counterpart of the branch of asymmetric 
solutions occurring in the case of two interacting, non-rotating Q-balls (boson stars). 

Decreasing 0J2, the numerical results show that the function fa progressively develops a local 
maximum on a ring in the equatorial plane. The function fa becomes progressively smaller, 
while the ring where the function fa reaches its maximum (typically \ fa\ ~ 1) becomes larger, 
spreading more and more over space. This phenomenon is illustrated in Fig|T7] (right) for 
Ui = 0.5, u>2 = 0.525 and 7 = —0.5. The pattern is qualitatively similar for other values of 
u)±, including uj\ = 0. Thus, the spherically symmetric, non-rotating Q-ball is progressively 
absorbed by the axially symmetric, rotating Q-ball when the parameter ui2 decreases. 
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4.1.2 Fixed lu 2 

Another striking feature of interacting Q-balls is that the direct interaction allows for solutions 
with arbitrarily small u% and hence charge Q\. Fixing L02 and varying u% shows that interacting 
solutions exist for lo\ > and that the function cf>i representing the non-rotating Q-ball spreads 
over the full r interval when oj\ tends to t-Oi yTnax . Only the rotating component of the solution 
(i.e. 4>2) survives in this limit. For 7 = —0.5 and UJ2 = 0.65 we find e.g. coi tmax « 0.75. 

4.2 Interacting boson stars 

It has been observed in [10] that for a single boson star the space-time becomes flat when the 
maximal value of oj is reached (consistent with the fact that the boson star vanishes and Q — > 
in this limit). In contrast, space-time become strongly curved close to the origin when the 
minimal value of u is approached. This feature is correlated to the fact that for small ui the 
matter fields become more concentrated in the region around the origin. 

We have studied the domain of existence for a = 0.1. This is given in Fig. [16] for 7 = 
(domain limited by lines with red bullets) and for 7 = — 1 (domain limited by lines with blue 
squares), respectively. The numerical integration becomes very difficult whenever approaching 
the limits of these domains of existence. The figure demonstrates that both gravity and the 
direct potential interaction leads to the existence of solutions with lower values of u\. We 
find that the interacting solutions stop to exist because the fields <j)\ vanishes identically when 
approaching a minimal value of L02 (and uj\ fixed). The other regions of the domain suggest the 
emergence of new branches of solutions which backbend from the branches we have constructed. 

To understand the features in more detail, we have studied the dependence of the properties 
of the solutions on the gravitational coupling a for 7 = 0, i.e. we let the boson stars interact 
via gravity only. We have chosen the frequency of the axially symmetric boson star to be fixed, 
L02 = 0.8, and have investigated the properties of the boson stars for varying a and frequency 
loi of the spherically shaped boson star. As discussed in more detail below, we observe that the 
number of solutions available depends strongly on u\ and that - in fact - additional branches 
of solutions exist. 

We have first studied solutions for which the frequency of the non-rotating boson star is 
smaller than that of the rotating boson star. We have chosen uj\ = 0.2, 0J2 = 0.8. For a = 0, 
no corresponding Q-ball solution exists. Increasing a we observe that a pair of solutions exist 
if a is sufficiently large. This is illustrated in Fig. [T8| where we give the mass M, the angular 
momentum J, /(0), <f>i(0) and <i>2,max in dependence on a. Clearly, two branches of solutions are 
present, which meet at a = a m i n ~ 0.3, such that for a < a m i n no solutions exist. The branch 
of solutions with lower mass likely exists for a — > 00, where the scalar fields slowly converge to 
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01 = 02 = 0. The branch with higher energy (we call it "the second branch") ceases to exist at 
a critical value of a = a cr , where /(0) is very small. Likely, further branches exist that start at 
the end point of this second branch. On these we would expect /(0) to decrease even further 
towards zero. 

1.6 -r 
1.4 - 
1.2 - 
1.0 - 
0.8 - 
0.6 - 
0.4 - 
0.2 - 
0.0 -■ 

0.2 0.4 0.6 0.8 1.0 

a 

Figure 18: The mass M, the angular momentum J, the values /(0) and <fti(0) as well as (j>2,max 
are given in dependence on a for a non-rotating boson star with u\ = 0.2 interacting with a 
rotating boson star with uj2 = 0.8. 

Next, we have studied the case where the frequency of the non-rotating boson star is equal 
to that of the rotating boson star, i.e. 002 = W2 = 0.8. Our results are shown in Fig. [19] In 
this case, Q-ball solutions exist in the flat space-time limit a = 0. For a > 0, a single branch of 
boson stars emerges from these flat space-time solutions. The mass M and angular momentum 
J of the interacting boson stars decrease with increasing gravitational coupling. The values of 
4>2,max and /(0) decrease, but stay finite for all values of the gravitational coupling. However, 
the value of 0i(O) tends to zero for a — > a cr ~ 0.24 indicating that cj)\ = 0. This means that for 
a > a cr the spherically symmetric boson star has disappeared from the system and the solution 
is simply a rotating, single boson star. 

Apparently, the pattern of solutions is quite different for the two cases studied. To understand 
the connection between the two patterns, we have studied the solutions for an intermediate value 
of uj\. We have chosen uj2 = 0.4 and 0J2 = 0.8. We observe that in this case, both patterns of 
solutions coexist. For a = both individual Q-balls - the non-rotating Q-ball and the rotating 
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Figure 19: The mass M, the angular momentum J, the values /(O) and (fii(0) as well as 4>2,max 
are given in dependence on a for a non-rotating boson star with uo\ = 0.8 interacting with a 
rotating boson star with uj2 = 0.8. Here 7 = 0. 

Q-ball - exist and do not interact (remember 7 = 0). Increasing a from zero, these solutions 
interact via gravity. As is shown in Fig. [20l the mass M and angular momentum J decrease for 
increasing a. For a ~ 0.110, however, W6 find, tllclt J — > cUld (j)2 max 

— ► 0. This is related to 

the fact that in this limit <p2 = such that the rotating boson star disappears from the system 
and the solution is simply a single non-rotating boson star. This is similar to the pattern found 
for lo\ = u)2 = 0.8 though it is the rotating boson star that disappears here. 

However, we find additional branches for sufficiently large a, which are disconnected from 
the a = solutions. These branches are shown in Fig. [2TJ One branch (the branch with lower 
energy) exists up to a — ► 00, while the second branch (the one with higher energy) emerges 
from the first branch at a > 0.113 and extends back in a. On this second branch, the value of 
/(0) decreases with increasing a such that at some critical value of a, /(0) becomes very small. 
Likely further branches exist on which /(0) tends to zero. This behaviour is similar to that 
observed for lui = 0.2 and u>2 = 0.8. 

To summarize, we find at least three branches of interacting solutions in this case: one 
branch that is directly connected to the flat space-time limit and exists for a S [0,0.110], a 
second branch that exists for a E [0.113, 00] and a third branch that exists for a £ [0.113, 0.185]. 
The second and third branch join at a = 0.113, while the first branch seems disconnect from 
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Figure 20: The mass M, the angular momentum J, the values /(O) and (fii(0) as well as 4>2,max 
are given in dependence on a for a non-rotating boson star with oo\ = 0.4 interacting with a 
rotating boson star with cc>2 = 0.8. Here 7 = 0. 

the others. 

Studying the equations for large values of the gravitational coupling a therefore suggests 
the following properties: (i) the interacting solutions obtained by deforming the flat space-time 
solutions via gravity do not stay bounded for large a and end up into single boson stars, where 
one of the two scalar fields vanishes identically, (ii) branches of solutions that have no flat space- 
time limit do not show the behaviour mentioned under (i) for large a. These branches can 
somehow be seen as non-perturbative since they do not exist for arbitrarily small values of a. 

5 Conclusions 

We have studied interacting Q-balls and boson stars, respectively, in great detail. While previ- 
ous work P2J 22j has focused on the different types of solutions possible, we have studied the 
parameter dependence of the mass and charges of the solutions here. 

We observe new features in the model. In the flat space-time limit the model describes 
interacting Q-haWs. We observe that when Q-balls are repelling, we can make the corresponding 
frequencies and charges arbitrarily small. This phenomenon appears for the choices of potential 
parameters done in previous work [5j [TO] . We would like to emphasize that single Q-balls with 
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Figure 21: The mass M, the angular momentum J, the values /(O) and <^>i(0) as well as 4>2,max 
are given in dependence on a for a non-rotating boson star with uo\ = 0.4 interacting with a 
rotating boson star with uj2 = 0.8. Here 7 = 0. 

arbitrarily small charges are also possible, but for different choices of the potential parameters 
than those done here and in [5l [10] . 

Whenever Q-balls and boson stars, respectively are interacting, a new type of solution is 
possible: one for which Q\ ^ Q2 while oj\ = 0J2- We interpret this as a symmetry breaking 
in the model. In the flat space-time limit, this is only possible when the Q-balls are repelling. 
Moreover, asymmetric solutions always have a companion solution in the form of a symmetric 
solution with Q\ = Q2- The asymmetric solution has much higher energy than the symmetric 
solution. In a concrete physical setting, we would thus expect the symmetric solution to be the 
generic solution. 

This changes, however, when boson stars are interacting solely via gravity. In this case, there 
is a frequency range in which only asymmetric solutions exist. This is an important observation 
when considering boson stars as dark matter candidates. 
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6 Appendix: The equations of motion 

Using suitable combinations of the Einstein equations we find the following expressions: 

„2« t , * * _ 1 1 f o„:„2 fl7 2 2 ■ ~~ - . -COS0 



r z d rr f + dttf = -7777 -2sin^ 9l z m 2 + ArfW r f + 2—fWgf + fd e ldgf 
2 j I \ sin & 

- 2l(r 2 (d r f? + (d e f?) + r 2 fd r ld r f + 4r sin 2 6l 2 md r m 

- 2 sin 2 9l 2 ((d r m) 2 + (d e m) 2 )) 

+ a {^T tt + fr 2 T rr + fT ee + g + T w - 2^T t ^ (27) 



r 2 d rr l + d u l 



1 1 



cos ( 



Urfld r l + A . 
2 j I \ sin 

+ a2l(r 2 T rr + T ed ) 



■fWgl - f(r 2 (d r l) 2 + {d e lf 



(28) 
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r 2 d rr m + dum = — --77 ( — 4/Zm + 4rfld r m + 6 C ° S ^ fldgm — 4l(r 2 d r fd r m 
2 j I \ s'mO 

+ dgfdgm) + 3f(r 2 d r ld r m + dgldgm) + Arlmd r f — 3rfmd r l) 

+ «^(mT w -rT ( ,) (29) 
sin ft 



- ^ (r 2 (d r ff + (a,/) 2 ) + (r 2 {d r l? + (del? 

+ f (rd r l + C £^d e l^j + - g (r\d r gf + (d e gf - rgd r g 

2a 2 

+ a {-2gr 2 T rr -2gTee + - M T -T w ) (30) 

sin 



r 2 d rr g + d ee g = sin 2 6[r 2 (d r mf + (d e m) 2 + m 2 - 2md r m 



where the non-vanishing components of the energy-momentum tensor read: 

I 2 
T tt = fV fa) - sin 2 6 >jm 2 V '(<£i,<fc) + w 2 2 + w 2 ^ 2 , - -m(fciu;i<ft + fou^) 

+ 9 1 2 + ^2) + - sin 2 9^m 3 (k 1 uj 1 (l) 2 1 + fow^) + sin 2 #4 m V?</>? + ^i) 

- ^(u, 2 fc 2 2 + u 2 k 2 cf> 2 ) + 1 sin 2 Qj 2 m\k 2 ^ 2 + fc 2 2 ) + ^((5 r 0i) 2 + (d r <p 2 ) 2 ) 

- sin 2 ^m 2 ((9 r 0i) 2 + (drfo) 2 ) + l^((a e 0i) 2 + (d e </> 2 ) 2 ) 

g r Ig 



1 • 2 Q m 
sm # — 

T l g 



-zn 2 O—Ufyfa) 2 + (d e <h?) (31) 



T rr = ~V{<t>u<h) + -H^fcp 2 + w 2 ^ 2 , + -%m(k 1 ui(j)l + k 2 u> 2 (pl) 9 r 2 + ^2) 



J_^ m 2 (fc 2 2 + ^ + { Q Ai) 2 + _ L^Q^f _ 1(^0 2 ) 2 



T ee = -r 2 ^y(</>i,0 2 ) +r 2 ^(c40f + w^ 

+ ^m 2 (fe 2 2 + ^ 2 )--4-(A: 2 2 + A: 2 V 2 ) 
/ sin # 

- r 2 {d r 4i? - r 2 {d r 42? + (d e <f>i) 2 + (^ 2 ) 2 (33) 



T re = 2(d r 0i)(d e 0i) + 2(a r ^ 2 )(a^ 2 ) 
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(34) 



T ttp = rsin 2 OjmV^,^) - r sin 2 B^m(J{§\ + uo\§\) - 2sin 2 d-^m 2 {k 1 ujx<\i{ + fou^l) 

+ ^-m{kl4>\ + k\<\>\) - ^ sin 2 6^m 3 (k 2 <pl + + 2/ciWi0 2 + 2fc 2 ^2</>2 

+ rsm 2 ^^^^) 2 + (d r </> 2 ) 2 ) + ^^((a^x) 2 + {defo?) (35) 

T w = -r 2 sin 2 (9-^y(0i,(/>2) + r 2 sin 2 fl^ (w 2 2 + w|0 2 ) + 2r sin 2 ^m^w^ 2 + fot^l) 



lm 2 (^ 2 + A: 2 2 )-r 2 sin 2 A 



+ fe 2 2 + fc 2 ^ 2 + sin 2 0-^m 2 (fc 2 </> 2 + fef^l) " r 2 sin 2 ^((<9 r 0i) 2 + (d r <p 2 ) 2 ) 
1 



- sin 2 0-((0^i) 2 + (0^) 2 ) (36) 

Finally, the Euler-Lagrange equations read (i = 1,2): 

11/ <9y A; 2 / 2 g 

r 2 d rr 4>i + = — -— I -r 2 l 2 g^— - 2 . i //g<fc + 2— (mfcj + nvi) 2 (pi 



+ Ir/Ro, • 2^/W„o, • /(r 2 ^/ (37) 
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